Elementary excitations of the symmetric spin-orbital model: The XY limit 
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The elementary excitations of the ID, symmetric, spin-orbital model are investigated by studying 
two anisotropic versions of the model, the pure XY and the dimerized XXZ case, with analytical and 
numerical methods. While they preserve the symmetry between spin and orbital degrees of freedom, 
these models allow for a simple and transparent picture of the low-lying excitations: In the pure XY 
case, a phase separation takes place between two phases with free-fermion like, gapless excitations, 
while in the dimerized case, the low-energy effective Hamiltonian reduces to the ID Ising model with 
gapped excitations. In both cases, all the elementary excitations involve simultaneous flips of the 
spin and orbital degrees of freedom, a clear indication of the breakdown of the traditional mean-field 
theory. 
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The impact of orbital degeneracy on the low-energy 
properties of Mott-Hubbard insulators is currently at- 
tracting a lot of attention following the progress in syn- 
thetizing and studying materials with these characteris- 
tics [0. It was already pointed out a long time ago by 
Kugel and Khomskii that such systems should have low- 
lying orbital excitations in addition to spin excitations 
More recently, it has been suggested that the inter- 
play between both degrees of freedom can have more dra- 
matic consequences. For instance, under suitable condi- 
tions the orbital degeneracy can enhance quantum fluctu- 
ations in the spin degrees of freedom and lead to gapped 
spin excitations even in the 3D case j3) . Another interest- 
ing situation is the SU(4) symmetric case ^||| where the 
system cannot choose locally between the configurations 
(spin singlet x orbital triplet) and (spin triplet x orbital 
singlet). Then the mean- field approach that decouples 
spin and orbital degrees of freedom on each bond can- 
not be a good starting point in that case since it violates 
basic SU(4) relationships between correlation functions 
on a given bond, as emphasized in Ref. ||. As a conse- 
quence, the traditional picture of relatively independent 
spin and orbital excitations must be abandoned. A clear 
picture of the low-lying excitations in such a case is still 
lacking though. 

In this Letter, we concentrate on the symmetric case. 
The basic model is the SU(4) symmetric Hamiltonian 
given by 



H = J 



2 Si • Si- 



2t, 



(1) 



where Si and are spin-1/2 operators corresponding to 
spin and orbital degrees of freedom respectively. This 
model has already been studied rather extensively by sev- 
eral methods [^-||. In particular, it is known from the 
Bethe ansatz solution that there are three branches of 
low-energy excitations . The physical interpretation of 



these branches is not straightforward though. The essen- 
tial complexity comes from the large local degeneracy: 
For a single bond, the groundstate is six-fold degenerate 
(spin singlet x any of the three orbital triplets or any of 
the spin triplets x orbital singlet). It is thus interesting 
to study the XXZ version of the model defined by 



x (2(rf rf +1 + t?t? +1 + Ar?r? +1 ) + -1 (2) 

In that case, the degeneracy is lifted within the triplet 
sector as soon as A < 1, and the groundstate of a given 
bond is only two-fold degenerate (spin singlet x orbital 
triplet with r t z ot = or spin triplet with S* ot = x orbital 
singlet). The essential ingredient, namely the symmetry 
between spin and orbital degrees of freedom, is preserved, 
but the Hilbert space of the low-lying sector is consider- 
ably reduced. 

In the following, we will concentrate on two versions 
of this model for which a transparent picture of the low- 
lying excitations can be obtained: 

The pure XY model: It corresponds to the previous 
Hamiltonian [Eq. (Q)] with A = and J; = J for all 
bonds, which can be written more compactly as 

H = JY, ( S t S i+i + S i S t+i) (rtrf+x + r^rt +l ) (3) 

i 

or in expanded form 

i 

1 + STr/te) ■ (4) 

i 

Analyzing this hamiltonian in the product basis ®i\rj)i, 
where 1 77) ^ = \Sf — ±i, r? = ±i) and denoting 
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|a±> 4 = ISf = ±1/2, 7f = ±1/2), 
|&±>i = |Sf = ±1/2, 7? =Tl/2), 



(5) 



one can easily show that all the matrix elements of H 
between states {|a±)j} and {|6±)j} vanish. As a conse- 
quence, the eigenstates can be classified according to the 
sequence of domains of phase A (with parallel orbital and 
spin states) and phase B (with antiparallel states) that 
include only states of type |o±)< and |6±)j respectively, 
and Eq. (Ji|) can be written as 



H = 2jJ2 + a % a+ +1 ) 

i 

+2J ^(PtPi+i + 



(0) 



where af — Sfrf" {(if = Sfrf) are the raising and low- 
ering spin-1/2 operators corresponding to |a±)j (|6±)j). 
Now the Hamiltonian within a phase can be mapped onto 
spinless fermions with a Jordan- Wigner transformation, 
and the groundstate energy of a domain of length L is 
given by 



E(L) = -2Jcos 



( (L/2+1)tt 
V 2(L + 1) 



sm 



2(L + 1) 



( 2(TTT) ) 



(7) 
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It is easy to check that E(L) < E{L{) + ... + E(L n ), 
where L\, L n are integers such that L% + ... + L n = L. 
So the model of Eq. (^|) undergoes a phase separation, 
and the groundstate is two-fold degenerate with a single 
domain of phase A or B. 

The phase separation has some drastic consequences 
for the thermodynamics, and the low-temperature prop- 
erties of the pure XY spin-orbital model turn out to be 
significantly different from those of the SU(4) symmetric 
model [Eq. (0)]. In Fig. ^ we show the entropies s per 
site as a function of temperature for both models, the 
SU(4) symmetric model [Eq. ([|)] and the pure XY-case 
[Eq. (0)]. They have been calculated numerically for 
chains of length L = 200 with periodic boundary condi- 
tions, using the continuous time quantum Monte Carlo 
loop algorithm At very low temperature, both 

entropies show a linear behavior (see inset of Fig. [I]), 
indicating the presence of gapless excitations. But the 
slope for the SU(4) is much larger than that in the pure 
XY spin-orbital case. In fact, in the first case the slope is 
three times bigger than that of single SU(2) antiferromag- 
netic Heisenberg chain j(| , while the slope of the entropy 
at low T in the pure XY spin-orbital model is equal to 
that of the XY model with coupling 2 J, as expected from 
Eq. (||). A further difference is visible as the tempera- 
ture is increased. The entropy of the SU(4) symmetric 
spin-orbital model remains approximately linear also in 
the intermediate temperature range (up to T w 0.2 J). 
The entropy of the pure XY spin-orbital model, on the 



other hand, coincides with the linear behavior of the XY- 
model only up to a crossover temperature T* m 0.05J. 
Above T* a sharp increase of the entropy takes place (see 
inset of Fig. |l|) . This increase comes from the additional 
entropy contribution of the domain walls which are more 
and more frequent with increasing temperature. This 
can be seen in Fig. ^ where the average density of do- 
main walls p (average number of domain walls per site) 
is depicted as a function of temperature. Up to T*, the 
spin-orbital model is in one of the phases A or B and 
p = within the statistical errors of the Monte Carlo 
simulation. Above T* the number of phase sectors in- 
creases very sharply with increasing T. 




FIG. 1. Temperature dependence of the entropy of the 
SU(4) symmetric spin orbital model of Eq. ([j]) (dashed line, 
taken from Ref. ||) and of the pure XY spin-orbital model 
of Eq. (^|) (solid line) for chains of length L = 200 with peri- 
odic boundary conditions in two different temperature scales. 
For comparison the analytical result for the entropy of the 
XY-model with coupling 2J and infinite length is also shown 
(dotted line). 

This behaviour can be easily understood in terms of 
the following, approximate free energy per site: 



f(T,p) = pE DW - Ts DW (T,p) - Ts D (T) 



(8) 



The various quantities entering this expression are: i) p, 
the concentration of domain walls, to be determined by 
minimizing the free energy; ii) Edw i the energy of a do- 
main wall. This is the energy required to split a finite 
chain of length L into two chains of length L — L\ and 
Li, i.e. E(L - Li) + E{Li) - E(L), where E{L) is de- 
fined in Eq. (Q). Edw is a priori a function of L and 
L\. It turns out that, for large enough L, Ejjw — 0.36 J 
regardless of L and L\ except for L\ < Lq ~ 20. Since 
Lq does not depend on L, this difference will play no role 
in the thermodynamic limit at low temperatures, and 
one can safely assume Edw = 0.36J; iii) sdw{T,p), the 
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entropy of the domain walls. For small p, it is given by 
Sdw{T,p) = —p\np; iv) sd(T), the entropy contribution 
of the domains. When p is small, finite-size effects are 
negligible, and sd(T) is equal to the entropy of the XY 
model with coupling 2J, i.e. (n/6)(T/J) at low temper- 
ature. 
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FIG. 2. Temperature dependence of the average density p 
of domain walls in the pure XY spin-orbital model. The error 
bars are smaller than the symbols. The inset shows Inp as 
function of the inverse temperature j3J as well as fits of the 
form p = exp(-S /T) and p = (L - 1) exp(-£i/T) in the 
intermediate and very low temperature ranges (for details see 
text). 

As for the lD-Ising model, minimizing with respect to 
the density of domain walls p leads to p — exp(— Eow/T) 
and s DW {T) = (E DW /T)exp(-E DW /T). The total en- 
tropy, which is the sum of sdw{T) and s_d(T), is then 
dominated at low temperature by sd(T) ~ (ir /6)(T/ J), 
while the domain wall contribution takes over at higher 
temperature. To be more quantitative, let us define the 
temperature T\ where both contributions are equal. It 
is given by (E DW /Ti) exp(-E DW /Ti) = (ir/6)(T 1 /J), 
which leads to T\ = 0.074J. This is in very good agree- 
ment with the numerical results of Fig. |l|. 

The prediction for the density of domain walls is also in 
very good agreement with the numerical simulations (see 
inset of Fig. ||, dashed line) for not too low low temper- 
atures: A fit with an exponential law p — exp(— Eq/T) 
for 8 < (3 < 12 gives Eq — 0.36(1), in very good agree- 
ment with the domain- wall energy Eow- For very low 
temperatures, namely for temperatures where the aver- 
age number of domain walls is of order 1 or smaller, and 
for finite systems, the above picture cannot work because 
the numerical simulations were performed using periodic 
boundary conditions, and domain walls can be created 
only by pairs with a minimum energy 2Eow- Neglect- 
ing configurations with more than one pair of domain 



walls, one can show that the concentration is expected 
to behave like p ~ (L - 1) exp(-2E DW /T). A fit of 
the very low temperature numerical data could indeed 
be performed with the law p = (L — 1) exp(— E\/T) for 
15 < < 25 (see inset of Fig. |]) with E x = 0.72(2), 
again in very good agreement with 2Eow- Note that 
the temperature To below which finite-size effects due 
to periodic boundary conditions start to influence the 
thermodynamics is given by exp(— Edw/Tq) = l/L, i.e. 
To ~ Ebw I hi L. Since it vanishes when the system size 
goes to infinity, the free energy of Eq. (f|) is expected to 
be valid down to zero temperature in the thermodynamic 
limit. 

To summarize this section, the very low temperature 
excitations correspond to simultaneous flips of spin and 
orbital degrees of freedom within one domain, and do- 
main wall excitations corresponding to collective exci- 
tations involving spins and orbital degrees of freedom 
play an important role above a cross-over temperature 
T* ~ 0.05J. 

The dimerized XXZ model: It corresponds to the 
Hamiltonian of Eq. @ with J. t — J if i is even and Ji = 
a J if i is odd. We wish to study that model in the limit 
a«l. Let us start by introducing some notations. The 
Hilbert space of a given bond is spanned by the 16 states 
\SS), {ST,} (t = -1,0,1), \TiS) (t= -1,0,1) and |T 2 T,) 
= —1,0,1), where the first (second) letter refers 
to the spins (orbitals), while \S) and |Tj) are the usual 
singlet and triplets given by \S) = (| Tl) - I it»/V2, 
\Ti) = | H>, |2b) = (I U) + I IT»/V2 and |T_ X ) = | ||>. 
In the pure Heisenberg case, the six states \STi) and \TiS) 
(i = —1,0, 1) are degenerate groundstates. However this 
degeneracy is partially lifted if A < 1 in Eq. (j2j) , and the 
groundstate of a given bond is only two-fold degenerate 
{\ST } and \T S}). If a = 0, the groundstate of Eq.(|) 
is then 2 L / 2 -fold degenerate, where L is the number of 
sites, since each dimer (i, i + 1), i even, can be in any 
of the two states \STo)i or \ToS)i. Let us study how 
this degeneracy is lifted when a is switched on. Since 
we have a two-level system on each dimer + 1), i 
even, we can define a pseudo spin-1/2 operator at that 
acts on this dimer with the identification \STq) = \ I) 
and \TqS) = \ f). An effective Hamiltonian can then be 
derived using standard many-body perturbation theory. 
The result depends on A. If A > 0, then the perturbation 
is lifted to first order in a, and the effective Hamiltonian 
reads 

H eff = «^ E (^+2 + \) ( 9 ) 
i even 

while if A = one has to go to second order perturbation 
theory to lift the degeneracy, and the effective Hamilto- 
nian reads: 

H e7f=~ a2 E« CT ?+2 + J) (10) 

i even 
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Several remarks can be made about these results: First 
of all, the effective Hamiltonian is always an Ising model 
to the first non-vanishing order in perturbation theory. 
Second, the effective Ising model is antiferromagnetic if 
A > and ferromagnetic if A = 0. These models are 
of course related by a simple transformation, but we ex- 
pect to have a transition line between these cases in the 
(a, A)-plane along which the effective Hamiltonian pre- 
sumably takes a more complicated form. Finally, and 
more importantly, we obtain an Ising model in terms of 
°f = i a t + cr^ ) / 2 , not erf. So the eigenstates must 
be written in terms of the eigenstates of erf, namely 
(\ST ) + \T S))/y/2 and (\ST ) - \T S))/y/2. They are 
thus of the general form 2'^^ even (|ST )i ± \T S)i). 

Let us now briefly discuss the low-energy properties. 
Quite generally, we expect to have a two-fold degenerate 
groundstate, and a gapped excitation spectrum. More 
specifically, the groundstates are given by 

\GS±) = 2- L / 4 J] (|5T )i ± \T Q S)i) (11) 

i even 

in the ferromagnetic case, and by 



\GS±) = 2~ L / i ][[ (\ST )i±(-lf 2 \T S) i 



(12) 



in the antiferromagnetic case. The first excitations are 
obtained by replacing (|ST ) - \T S))/V2 (resp. (\ST ) + 
\T S))/V2) by (\ST ) + \T S))/V2 (resp. (\ST ) - 
\TqS))/ v2) in one of these groundstates with energy a 2 /2 
in the ferromagnetic case and a A 2 /2 in the antiferromag- 
netic case. So there is indeed a gap in the spectrum. 
More importantly, it is clearly impossible to separate spin 
and orbital degrees of freedom since (|ST ) + \T a S))/\/2 
and (\ST )~\T S))/V2a,re not eigenstates of (Si+S i+ i) 2 
or of {fi +fi+i) 2 , and the excitations are neither spin ex- 
citations nor orbital excitations. They are transitions 
between resonating valence-bond states that intimately 
mix spin and orbital degrees of freedom. 

It is also interesting to note that the correlation func- 
tions on a strong bond (i even) (Si.Si+%), (7^.7^+1) and 
((Si.Si + i)(fi.T i+ i)) are all negative, as in the SU(4) sym- 
metric case, which excludes mean-field theory as a good 
starting point for the same reasons (see Ref. J(|). 

Discussion: Coming back to the original problem of 
the nature of the excitations in the SU(4) symmetric 
model, let us put our results in perspective. In both cases 
studied above, exact results have been obtained, and the 
low-energy excitations are neither spin nor orbital ex- 
citations, but involve both spin and orbital degrees of 
freedom on an equal footing. This is a clear indication of 
the breakdown of mean-field theory. It strongly suggests 
that the model of Eq. (Q) also possesses such low-lying 
excitations. In particular, the operators of = S^t^ 1 and 
case are linear combinations of 
the operators S a T^ which can be seen as generators of 



the SU(4) algebra (see Ref. ||), and it is likely that at 
least part of the low-lying modes of the SU(4) symmetric 
model will be predominantly built on these generators 
and will retain the mixed character observed here. Be- 
sides, the fact that the correlation functions (Si-Si+i), 
{n-n+i) and ((Si.Si + i)(n.T l+1 )) are all negative on a 
given bond appears in the XXZ case as a direct conse- 
quence of the local degeneracy between the states (spin 
singlet x orbital triplet) and (spin triplet x orbital sin- 
glet). So the picture that emerges is that the symmetry 
between spin and orbital degrees of freedom has dramatic 
consequences on the low-lying excitations: The system 
is not able to choose between spin or orbital singlets or 
triplets, and the excitations are an intricate mixture of 
spin and orbital degrees of freedom. 

To complete the picture, it will be useful to study the 
XXZ model in the whole parameter range < a, A < 1. 
The main issues are: i) The evolution of the spectrum 
along the line (a = 1, A = 0) — > (a = 1, A = 1) joining 
the XY case and the model of Eq. (|l|) ; ii) The number of 
low-lying modes, and in particular the presence of a gap 
as a function of a and A; iii) The nature of the effective 
model in the limit A ~ 1, a <C 1. Work is in progress 
along these lines. 
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